Topological px + ipy Superfluid Phase of a Dipolar Fermi Gas in a 2D Optical Lattice 
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In a dipolar Fermi gas, the anisotropic interaction between electric dipoles can be turned into 
an elfectively attractive interaction in the presence of a rotating electric field. We show that the 
topological Px+ipy superfluid phase can be realized in a single-component dipolar Fermi gas trapped 
in a 2D square optical lattice with this attractive interaction at low temperatures. The px + ipy 
superfluid state has potential applications for topological quantum computing. We obtain the phase 
diagram of this system at zero temperature. In the weak-coupling limit, the p-wave superfluid phase 
is stable for all filling factors. As the interaction strength increases, it is stable close to filling factors 
n = or n = 1, and phase separation takes place in between. When the interaction strength is 
above a threshold, the system is phase separated for any < n < 1. The transition temperature of 
the Px + ipy superfluid state is estimated and the implication for experiments is discussed. 
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Unconventional quantum states of ultra-cold Fermi 
gases have attracted a lot of attention in the recent theo- 
retical and experimental studies. Among them, one of the 
most desirable states is the topological superfluid Px+iPy 
phase in a two-dimensional (2D) single-component fermi 
gas. The zero-energy Majorana modes were predicted 
to appear inside vertex cores in this p-wave state 
Due to their insensibility to local perturbations, Ma- 
jorana fermions are very appealing for applications in 
topologically-protected quantum computation [2] . 

The p-wave superfluid phase was predicted to appear 
near a p-wave Feshback resonance in a Fermi gas [3, 3] . 
However the lifetime of the Fermi gas is severely lim- 
ited by three-body collisions, and so far this superfluid 
regime has not been achieved in experiments. The suc- 
cessful creation of '*''K^''Rb polar molecules ^] has pro- 
vided another platform to study unconventional quantum 
states of Fermi systems. For dipoles aligned parallel to 
the z-axis, a p-wave superfluid state with the dominant 
Pz symmetry was predicted in a three-dimensional dipo- 
lar fermi gas Q . In a dipolar fermi gas trapped in x — y 
plane, a Px or py superfluid state may appear if the angle 
between the dipole direction and the x — y plane is small 
enough [t^ . The Px + ipy superfluid phase was predicted 
for a 2D dipolar Fermi gas with a circularly-polarized 
microwave field near the resonant frequency (8|. The d- 
wave and extended s-wave superfluid phases were also 
predicted for a similar setup in an optical lattice with 
a linearly-polarized resonant microwave field [9]. In ex- 
periments, however, the precise control of the microwave 
polarization has been difficult ^lO^. 

Here we propose a different method to realize a stable 
Px + iPy topological superfluid phase in a dipolar Fermi 
gas. For *°K^^Rb or other fermionic polar molecules 
trapped in a 2D square optical lattice, the dipole direc- 
tion can be fixed by applying an electric field. If the 
electric field makes a small angle with the lattice plane 
and rotates fast around the axis normal to plane, the 
time-averaged interaction between dipoles is effectively 



attractive |ll| . The px + ipy topological superfluid phase 
can be realized in this dipolar Fermi gas at low temper- 
atures due to the attractive interaction. 

This work is organized as follows. First we derive 
the effective Hamiltonian for a single-component dipolar 
Fermi gas trapped in a square lattice with a fast-rotating 
electric field. Then we study the Px +ipy superfluid state 
in the self-consistent Hartree-Fock approximation. The 
phase diagram of this system is obtained at zero tempera- 
ture. When the interaction is weak, the p-wave superfluid 
state is stable near filling factors n — and n — I, and 
phase separation occurs in between close to half filling. 
When the interaction strength is above a critical value, 
the system is phase separated for any < n < 1 . We also 
discuss the superfluid transition temperature and its im- 
plication for experiments in the end. 

Effective Model. We consider a single-component dipo- 
lar Fermi gas trapped in a square optical lattice under a 
rotating electric field 

E(f) — E[cos (pz + sin (p{cos fltx + sin flty)] , (1) 

where fl is the rotation frequency, E is the magnitude of 
the electric field, and (p is the angle between the electric 
field and z direction. The optical lattice is in the x — 
y plane and the rotation axis is along the z direction. 
In strong electric fields, dipoles are aligned parallel to 
E(i). For fast rotations, the effective interaction between 
dipoles is the time-averaged interaction 
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where d is the dipole moment, and r is the distance be- 
tween the two dipoles. In the following we consider the 
case with cos ip < where the effective interaction is 

attractive, Vdd{r) < 0. 

The effective Hamiltonian of this dipolar Fermi gas is 
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given by 

<ij> i-^j 

where Ci is the fermion annihilation operator at lattice 
site i — (ixyiy)-, Vi-j is the interaction between fermions 
at site i and j given by Vi-j = Vdd{\^i — and is 
the lattice vector. For simplicity we assume that in the 
optical lattice the fermion can only hop between near- 
est neighbors with hopping amplitude t. The interaction 
strength is the largest for nearest neighbors, denoted by 
J = |Vdd(a)| = d^|3cos'^(y5 — l\/{2a^), where a is the lat- 
tice constant. 

The Hamiltonian in Eq. ([3]) has the important particle- 
hole symmetry. Under the transformation c| Ci and 
Ci — !> c|, a state with filling factor n turns into another 
state with filling factor 1 — n. The chemical potential fi 
transforms as /i — > — /i + V'(O) where V{0) = J2i^o ^Oi — 
—9.02 J. Note that the hopping amplitude t becomes —t 
under this transformation, but in the square lattice this 
phase change can be eliminated by a unitary transfor- 
mation Ci — —Ci for only one sublattice with the same 
odevity of ix + iy 

I 
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where = Eh + + Ep. 
In momentum space, we have 

H' - ^iNf = 5]Kk4ck + ^c_kCk + ^clc\] ~ ^c, 

k 

(5) 

where ■fk = Ek + — ^, £k = —2t(coskxa + cos kyu) is 
the band energy, /i is the chemical potential, Sk is the 
Hartree-Fock self-energy given by 

Sk- V(0)n-l^y(k-k')nk', 

k' 

the pairing gap Ak is given by 

k' 

N is the total number of lattice sites and Nf is the 
fermion number operator, 

y(k) = ^Vm exp(-ik • r„j), 
5k = (c-kCk), and ny, = (cJ^Ck). 



Self- consistent Hartree-Fock approximation. Due to 
the attractive interaction, fermions tend to pair with each 
other at low temperatures. In the square lattice, since 
the interaction strength is the largest between nearest 
neighbors, and the pairing probability is the largest for 
fermions from nearest neighbors. Under the constraint of 
Fermi statistics, the dominant pairing symmetry should 
be p-wave. To study this superfluid pairing state, we use 
the self-consistent Hartree-Fock approximation in which 
the interaction produces the Hatree energy Eh , Fock en- 
ergy Ex, and the pairing energy Ep given by 

E- = -\T.^^~3(Ac3){c]c,), 

Ep - \ Y,V,^,{c\c]){c,c,), (4) 

where n = (c|ci) is the filling factor. The Hamiltonian 
in Eq. ^ can be approximated as 



•]) - {{4c3)c]cz + h.c.) + ((cJc])cjQ -t- h.c.)] - Ec, 

I 

The Hamiltonian in Eq. ([5]) can be diagonalized by Bo- 
goliubov transformation. The Hartree-Fock self-energy 
and gap can be determined self-consistently 

Ek = F(0)n-l^F(k-k')i[l-^tanh(^i?k')], 

k' 

(6) 

= -^E^(k-k')^tanh(^i?k'), (7) 

where the quasi-particle excitation energy E^ is given by 
El, = |Akp and /3 = l/iksT). The filling factor 

n can be also determined self-consistently. 




We numerically solve the Hartree-Fock self-energy 
equation ([6]), gap equation ([7]), and number equation (|8]) 
together, and find that the gap Ak has finite magnitude 
and displays a complex symmetry below a critical tem- 
perature. A global unitary transformation can always be 
applied so that the real part of the order parameter Ak 
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FIG. 1: The gap Ak versus wavevector k at J/t = 0.6, 
n — 0.11, and T = 0. The solid and dashed lines are the 
real and imaginary parts of the gap. (a) At kya/n — 0, the 
real part of gap has the same sign structure as the function 
sm{kxa), and the imaginary part vanishes, (b) At kxa/n = 0, 
the imaginary part of gap has the same sign structure as the 
function sin(fcya), and the real part vanishes. 



remains the same under the transformation ky — —ky 
but changes sign when kx —kx, while the imaginary 
part of the gap transforms oppositely. As shown in FiglTJ 
we find from numerical computation that the real part 
of the gap has the same sign structure as the function 
sm{kxa) and the imaginary part has the same sign struc- 
ture as the function sin(/cya). 

The magnitude of the gap is a monotonically increas- 
ing function of the interaction strength J. For the fixed 
interaction strength, the magnitude of the gap also in- 
creases with the filling factor up to half filling, as shown 
in Fig. [2] Due to the particle-hole symmetry, one state 
at filling n — 5 is SL dual to another state at n — 1 ~ S, 
and they share the same order parameter. Thus the gap 
of the superfluid state above half filling can always be 
obtained from the state below half filling. At half-filling, 
the magnitude of the gap is the largest. 

Zero- Temperature Phase Diagram. Due to the attrac- 
tive nature of the effective dipole-dipole interaction, the 
free energy of this dipolar Fermi gas is smaller than that 
of an ideal Fermi gas trapped in the same optical lattice. 
This energy reduction caused by the interaction increases 
with the filling factor. When the interaction strength is 
large enough, the interaction effect is dominant and the 
system can be unstable. As shown in Fig. [3l at T = 
and J/t = 0.7, the chemical potential is a monotonically 
decreasing function of the filling factor and the compress- 
ibility is negative for 0.3 < n < 0.7, indicating that the 
Px + ipy superfluid state is dynamically unstable. The 
system is phase separated in the unstable region and be- 
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FIG. 2: The gap magnitude Ao versus filhng factor n for 
J/t = 0.6 and T = 0, where Aq = |Ak| at k = {n/2a,n/2a). 
At half fiUing, Ao has the largest value. The gaps for fiU- 
ing factors 0.5 < n < 1 can be obtained from particle-hole 
symmetry. 




FIG. 3: Chemical potential fi versus the filling factor n at 
J/t — 0.7 and T — 0. The dashed line shows the phase- 
separation region. The shaded areas above and below the 
dashed line are equal, as required by the phase-separation 
condition. 



comes a mixture of n = 0.12 and n ~ 0.88 states, shown 
by the dashed line in Fig. [3l Due to the particle-hole 
symmetry, the phase-separation region is symmetric with 
respect to the half-filling point. 

We find that the phase-separation region increases with 
the interaction strength. In the weak-coupling limit, 
J -> 0, the phase-separation region vanishes. With 
strong couplings, for any < n < 1, the px + ipy super- 
fluid state is unstable and the system becomes a mixture 
of n = and n = 1 states which are insulating states due 
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FIG. 4: Phase diagram at zero temperature. The sohd hne is 
the boundary between the p^c + ipy superfiuid phase (p-wave 
SF) and phase-separation region (PS). The superfiuid state 
vanishes when J > 0.89t 



to zero particle or hole densities. At the critical inter- 
action strength where the superfiuid state vanishes, the 
chemical potential satisfies /i(0) = where fi{n) is 

the chemical potential at filling factor n and /i(0) — —At. 
Due to the particle-hole symmetry, /x(0) + ^(1) — V{Q), 
the critical interaction strength satisfies V{Q) = —8t or 
J — 0.89t. When the interaction strength is above the 
critical value J = 0.89^, the system becomes a mixture 
of n = and n = 1 states. The zero-temperature phase 
diagram of this system is obtained based on these results 
and shown in Fig. H) 

Discussion and Conclusion. When the temperature 
increases, the phase-separation region decreases due to 
the increase in entropy. For example, we find that at 
ksT = O.lt and J = t the phase separation nearly disap- 
pears and the superfiuid phase exists near half filling. In 
the superfiuid region, as the temperature increases, this 
2D dipolar Fermi gas eventually undergoes a Berezinskii- 
Kosterlitz-Thouless (BKT) transition from the superfiuid 
state to the normal state. At the critical point, the 
vortex-antivortex pair disassociates and it costs zero free 
energy to generate a single unbound vortex. 

The transition temperature Tbkt is given by [l5| , 

Tbkt = ^P, (9) 

where p is the superfiuid density. In the Px + ipy super- 
fiuid state, the superfiuid density refers to the geometric 
mean of eigenvalues of the superfiuid density tensor which 
describes the response of the system to phase twists in 
the superfiuid order parameter. At finite temperatures, 
the superfiuid density is reduced due to excitations of 
fermionic quasi-particles and bosonic (collective) modes. 
Here, we assume that the bosonic excitations are negli- 



gible and the superfiuid density tensor is thus given by 
P^3 = ^{nk^^^j£_i, - lka,^k9j^k) (10) 

k 

where ^^f = df/dh and Fk = /3sech2(/3£;k/2)/4 is 
the Yoshida distribution. Since the off-diagonal ma- 
trix elements of the superfiuid density tensor are small, 
the superfiuid density is approximately given by p w 
(Pxx + Pyy)/'^- By solving Eq. (jO]) numerically, we obtain 
the BKT temperature ksTBKT = O.llt for J = t and 
n = 0.48. 

In experiments, the hopping amplitude is approxi- 
mately t = lOOTi Hz 13] for a typical optical lattice with 
a = 1/zm. The strength of the effective dipole-dipole in- 
teraction can be tuned up to J = 200?i Hz for singlet 
rovibrational state of KRb molecules with experimental 
conditions 0, [l2| . The potential strength and lattice con- 
stant of the optical lattice can be tuned to increase or de- 
crease the hopping amplitude t and interaction strength 
J, which can be helpful for observing the Px + ipy super- 
fluid state. 

In summary, we propose that the topological Px + ipy 
superfiuid state can be created in a dipolar Fermi gas 
trapped in a 2D optical lattice in the presence of a rotat- 
ing electric field, due to the effectively attractive dipole- 
dipole interaction. The phase diagram of this system is 
obtained at zero temperature. The px + ipy superfiuid 
state is stable for any filling factor in the weak-coupling 
limit. As the interaction strength increases, phase sep- 
aration occurs near half filling. When the interaction 
strength is above a threshold, the superfiuid state is un- 
stable for any filling factor < n < 1 and the system 
becomes a mixture of n = and n = 1 states. We esti- 
mate the BKT transition temperature of the topological 
Px-'ripy superfiuid state which may be accessible in future 
experiments. We would like to thank T.-L. Ho for help- 
ful discussions. This work is supported by NSFC under 
Grant No 10974004. 
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